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impedance tomography system comprises a plurality 
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electrical conductivity distribution of which is to be 
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each pair of electrodes will be stimulated in turn by 
applying a known current to them. For each stimulated 
pair, a resulting potential difference will be recorded 
at the pair of the remaining electrodes which, if 
stimulated, would produce an electric field with vectors 
of the electric field produced by the stimulated pair of 
electrodes. 
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ELECTRICAL IMPEDANCE TOMOGRAPHY METHOD AND ELECTRODE 



ARRANGEMENT FOR USE THEREIN 



DESCRIPTION 

5 TECHNICAL FIELD: 

The invention relates to a method of performing electrical impedance tomography 
measurements and electrode arrangements for use therein. 

The invention is applicable generally to the imaging of electrical impedance 
variations in conductive media, and especially in volume conductors such as parts of the 
10 human body, the ground, pipelines, and so on. 



BACKGROUND ART: 

Electrical impedance tomography (EIT), also known as "applied potential 
tomography" (APT), is used to provide images of spatial variations in electrical 

15 impedance within a conductive volume conductor. The electrical impedance changes 
are measured by providing an array of electrodes about the volume to be imaged. A 
stimulus current is supplied to each pair of electrodes in turn and the resulting potential 
differences recorded between pairs of the remaining electrodes. The process is repeated 
until all of the independent combinations of stimulation/recording are exhausted. The 

20 measurements are used to determine the transfer impedance changes and construct an 
image of the volume. In contrast to X-ray computer tomography applications, where the 
paths of photons through a body are straight lines, the current paths in EIT are functions 
of an unknown electrical conductivity distribution. This gives rise to a non-linear image 
reconstruction problem. Nevertheless, algorithms have been devised for converting the 

25 series of measurements into an image of the electrical impedance distribution within the 
body at a sufficiently high rate that changes associated with respiratory and cardiac 
functions, for example, can be monitored. 

The quality of the image depends upon the sensitivity of the measuring system, 
which typically varies with distance from the electrodes. In clinical applications, for 

30 example, where the electrodes are distributed around the body part to be measured, such 
as the thorax, difficulties may be encountered in resolving impedance changes which 
occur deep in the body. Most attempts to improve image quality have focused upon the 
reconstruction algorithms, as disclosed, for example, in U.S. patents Nos. 4,617,939, 
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5,381,333 and 5,465,730. Limited attention has been devoted to the electrode 
arrangement. For example, a study of various electrode configurations was described 
by Booth et ah in "A Comparison of Three Electrode Configurations for Electrical 
Impedance Tomography", IEE Colloquium on "Innovation in Instrumentation for 
5 Electrical Tomography", pp. 11/1, 11/3, 1995. British patent application number 
2,257,530 disclosed electrodes in a circular array or "rosette" placed upon a substantially 
flat part of the body, such as the chest. This was said to be especially desirable when 
imaging the heart, since impedance changes occurring in the heart during the cardiac 
cycle will be larger than those occurring in the region of the skin. 

10 Little attention has been paid to the way in which the measurements themselves 

are made. There are two kinds of procedure for taking the measurements. One involves 
applying stimulation currents simultaneously to a ring of electrodes to generate a current 
pattern and simultaneously measuring the corresponding voltage distribution around the 
same or adjacent electrodes. An example of such an approach is described in an article 

15 by P. Hua et ah entitled "An Electrical Impedance Tomograph using Compound 
Electrodes", presented at the IEEE Engineering in Medicine & Biology Society 11th. 
Annual International Conference, 1989. Hua et ah described experiments with thirty-two 
compound electrodes encircling a body to be measured. Each compound electrode 
comprised an inner electrode surrounded by an outer electrode with an annular space 

20 between them. In some cases, Hua et ah short-circuited the inner and outer electrodes 
and used them for both stimulation and recording. In other cases, Hua et ah injected 
a spatially-sinusoidal current pattern into the ring of 32 outer electrodes and recorded the 
corresponding generally sinusoidal voltage distribution around the ring of 32 inner 
electrodes. This approach requires 32 adjustable current generators which must be 

25 adjusted individually to give the required pattern as accurately as possible, making it 
difficult to obtain repeatable measurements accurately. Also, the equipment is complex 
and costly. 

An example of the other kind of procedure is described in US patent number 
4,617,939 (Brown et ah) issued October 1986. Brown et ah describe positioning 16 
30 electrodes around a body, applying a current to a first pair of the electrodes and 
recording the potential difference between every other pair of the remaining electrodes. 
As also noted by Hua et ah {supra), significant contact impedance between the electrodes 
and the body surface militates against measurement sensitivity, so Brown et ah avoid 
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using the same electrode for both stimulation and recording. Brown et aL repeat the 
procedure, applying current to each pair of electrodes in turn. For each pair, Brown et 
al. measure voltages at all of the remaining electrodes. In practice, such a procedure 
would not be entirely satisfactory because many of the voltage measurements would be 
5 comparable with noise levels. 

The present invention seeks to mitigate the disadvantages of these known EIT 
systems and provide an improved electrical impedance tomography procedure with 
enhanced sensitivity to changes in the electrical conductivity distribution inside the body 
of interest. 

10 

DISCLOSURE OF INVENTION: 

According to one aspect of the invention, a method of determining electrical 
impedance tomography of a conductive volume comprises the steps of selecting, in turn, 
a plurality of pairs of locations upon a surface of the conductive volume, supplying a 

15 stimulation current to the surface by way of the selected pair of locations, and recording, 
for each selected pair of locations, a resulting potential difference between at least one 
pair of the remaining locations, characterized in that , for a particular pair of locations 
stimulated, the potential difference is recorded for a number of pairs of said remaining 
locations that is less than the total of said plurality of locations and also comprises that 

20 pair of locations which, if stimulated, would produce an electric field with vectors most 
closely aligned with the corresponding vectors of the electric field produced by 
stimulation of said particular pair of locations. 

The above method may be performed using a plurality of electrode means, each 
disposed upon or adjacent said surface at a respective one of said plurality of locations, 

25 each electrode means being used to apply stimulation current to, or record potential at, 
said respective one of said locations. 

Each electrode means may comprise two electrodes, that are closely located 
spatially, one for stimulation and the other for recording. One electrode may be 
disposed inside the other in a plane which, in use, will be parallel to the said surface. 

30 The stimulation current could then be applied to two outer electrodes and the potential 
difference measured at a said selected pair comprising the corresponding two inner 
electrodes. Conversely, the stimulation current could be applied to the two inner 
electrodes and the potential difference measured at the corresponding outer electrodes. 
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Alternatively, each electrode means may comprise a first electrode and a second 
electrode spaced apart from each other in a direction which, in use, will be substantially 
normal to the said surface. The stimulation current could be applied by way of two of 
said second electrodes and the potential difference measured between two of said first 
5 electrodes, or vice versa. 

According to another embodiment of the invention, a method of measuring 
electrical conductivity distribution within a conductive volume includes the steps of 
positioning adjacent a surface of the body an array of electrodes, with a first group of 
the electrodes closer to the surface than a second group of the electrodes, applying a 

10 stimulation current to selected pairs of electrodes in turn and measuring, for each pair 
stimulated, corresponding electrical potential differences produced between pairs of the 
remaining electrodes, and processing the measured potential differences to determine 
electrical conductivity variations within the body. 

According to another aspect of the present invention, an electrode arrangement 

15 for an electrical impedance tomography system comprises a plurality of electrodes and 
support means for supporting the electrodes adjacent a surface of a conductive body the 
conductivity of which is to be mapped by the system, the arrangement being such that 
a first group of said plurality of electrodes will be closer to the surface than a second 
group of the plurality of electrodes. 

20 According to yet another aspect of the invention, an electrode arrangement for 

an electrical impedance tomography system comprises a plurality of electrodes in an 
array mounted in a support medium for supporting the electrodes adjacent a surface of 
a volume to be imaged, the support medium having anisotropic conductivity, its 
conductivity in a direction that, in use, is normal to the surface, being significantly 

25 greater than its conductivity in a transverse direction. 

According to a further embodiment of the invention, an electrode arrangement for 
an electrical impedance tomography system comprises an array of electrodes fixed 
spatially relative to each other in a rigid support medium, and an interface medium for 
interfacing the rigid support to a surface of a conductive volume to be measured, the 

30 interface medium being conductive so as to connect the electrodes electrically to the 
surface and pliable so as to conform to variations in relief of the surface. 

An advantage of this electrode arrangement of such further aspect, as compared 
with electrodes which are simply attached directly to the volume, is that the relative 
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position of the electrodes, and the dimensions of the array, are fixed, which facilitates 
static image reconstruction of the image. 

According to yet a further aspect of the invention, an electrode arrangement for 
use in electrical impedance tomography comprises a plurality of electrode means each 
5 formed by a stimulation electrode and a recording electrode closely located spatially. 
In one embodiment of this aspect, the electrode means comprises an inner electrode 
surrounded by an outer electrode with an annular channel between them. Electrically 
insulating sealing means may be provided in the annular channel for contacting said 
surface, in use, and insulating the inner electrode electrically from the outer electrode. 
10 Various objects, features, aspects and advantages of the present invention will 

become more apparent from the following detailed description of preferred embodiments 
of the invention, taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS: 
15 Figure 1, labelled PRIOR ART, is a diagram of the generic components that are 

part of a typical EIT apparatus; 

Figure 2 is a diagram of a first embodiment of the invention comprising an 
annular EIT electrode array; 

Figure 3 is a three-dimensional view of the annular electrode array of Figure 2; 
20 Figure 4 is a cross- sectional detail view taken on the line IV-IV of Figure 2; 

Figure 5 compares the sensitivity values of the annular EIT electrode array of 
Figure 2 with those of the EIT system of Figure 1, for radial positions; 

Figure 6 illustrates the sensitivity theorem used in determining electrical 
impedance tomography (EIT) of a volume conductor; 
25 Figure 7 is a flowchart representing a typical measurement procedure using the 

electrode means of Figure 2; 

Figure 8 is a plan view of an alternative form of electrode means, comprising an 
inner electrode and an outer electrode (also referred to as a compound electrode; 
Figure 9 is a cross- sectional view of the electrode means of Figure 8; 
30 Figure 10 illustrates a second embodiment of the invention in the form of a planar 

rectangular electrode array for a ground probing application; 

Figure 11 illustrates one cell of the electrode array of Figure 10; 



WO 98/23204 PCT/CA97/00876 

6 

Figure 12 illustrates a third embodiment in the form of a cylindrical electrode 
array for three-dimensional imaging of a human torso; and 

Figures 13 A to 13J illustrate the image reconstruction obtained from a 
conventional EIT electrode array, such as that shown in Figure 1, and an electrode array 
5 in which the stimulating electrodes and recording electrodes are co-located, such as that 
shown in Figure 2. 



BEST MODES FOR CARRYING OUT THE INVENTION: 

Referring to Figure 1, a typical known EIT apparatus comprises an electrode 

10 array comprising sixteen electrodes B x to E 16 distributed, equally spaced apart, around 
a conductive volume 12, to which they are attached directly, if necessary with a suitable 
conductive gel ensuring good electrical contact. The conductive volume 12 may be, for 
example, a human thorax. The electrodes are connected by a link 14 to electronic 
conditioning circuitry 16 which will usually include components such as a current source 

15 controlled by a D-to-A-converter, an A-to-D converter, differential amplifier, filter, 
analog multiplexer, clock, and the digital I/O required to control the circuitry from a 
computer 18 to which the electronic conditioning circuitry is connected. The computer 
18 may be a personal computer equipped with a digital signal processor card used for 
the image reconstruction process and a suitable display 20 for displaying the images. 

20 Using the procedure disclosed in US 4,617,939, a first pair of the electrodes Ej 

and E 2 are stimulated by injecting a current of, say, 4mA at 50 kHz, into them and the 
corresponding potential differences measured between each of the other pairs of 
electrodes, i.e. (E^), (E 4 -E 5 ), (E 5 -E 6 ), and so on up to (E 15 -E 16 ). The potential 
differences between electrodes E 16 and E l5 B x and E^ and E^, and E3 are not measured 

25 because the voltage drop across the contact resistance due to the injected current would 
affect the measurements. 

The next pair of electrodes E 2 and E3 then are stimulated and the potential 
differences between pairs of the other electrodes measured. The procedure is repeated 
until all of the pairs of electrodes have been stimulated in turn. The accumulated 

30 measurements are processed by the computer 18 and an image of the electrical 
impedance variation across the conductive volume is created and displayed on display 
unit 20. Various algorithms have been disclosed for computing these impedance 
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variation, such as the so-called "back-projection" reconstruction algorithm and the 
regularized inverse of the sensitivity matrix. 

While this procedure might appear satisfactory theoretically, in practice, noise 
may obscure many of the readings and reduce the quality of the reconstruction. 

5 Embodiments of the present invention favour high signal-to-noise ratio measurements by 
limiting the measurements, for a given stimulation pair, to that pair or pairs whose 
"recording" electric field vectors most closely align with the electric field vectors of the 
stimulated pair of electrodes. It will be appreciated that the recording electrodes measure 
potential difference so their "recording" electric field vectors is really that which would 

10 be produced if a stimulation current were injected into them. Generally, the present 
invention aims to keep the stimulating and recording electric field vectors as close to 
each other as possible. (Ideally, they would be superimposed). Various arrangements 
for doing so will now be described. 

A first embodiment of the invention will now be described with reference to 

15 Figures 2 and 3 in which the electronic conditioning circuitry 16', computer 18' and 
display unit 20' are similar to those of Figure 1 and so have the same reference 
numerals, but with a prime. 

The electrode array 22 shown in Figures 2 and 3 differs from that shown in 
Figure 1 in that it comprises two groups of electrodes arranged in an annular 

20 configuration, i.e. in two rings, one within the other, and supported by an annular 
support. The inner group comprises recording electrodes Rj to R 16 and the outer group 
comprises stimulating electrodes Sj to S 16 . The annular support comprises annular 
flanges 24 and 26 (see Figure 4) bonded to opposite sides of an array of radial separators 
28j - 28 16 and a circumferential band 30 is bonded to their respective outer edges, 

25 forming a ring of wedge-shaped compartments or cells 32! to 32 16 which are open at their 
inner ends 34j to 34 16 and alternate with separators 28! to 28 i6 . The annular flanges 24 
and 26, separators 28, and outer band 30 are of insulating material. Each pair of an 
inner electrode R and an outer electrode S mounted in a corresponding one of the 
plurality of wedge-shaped compartments or cells 32j to 32 !6 . 

30 The wedge-shaped compartments 32 1 to 32 16 are equally distributed and each 

subtends a 15 degree angle. They contain a conductive material whose electrical 
conductivity is comparable to the expected average electrical conductivity of the imaging 
medium. For example, for measuring a human body part with a conductivity of, say, 
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0.002 mho/cm, the compartments 32 t to 32 16 might contain agar jelly formed by mixing 
agar powder with saline solution to form a conductive gel, giving a conductivity for the 
wedge-shaped compartment of about 0.002 mho/cm. The separators 28! to 28 16 are 
poorly conductive and each subtends an angle of 7.5 degrees so that the annular support 
5 is circular. The separators 28j to 28 16 might be of synthetic plastics material, for 
example polypropylene. 

In use, a conductive interfacing medium 36, for example sponge material soaked 
in saline solution, is provided between the conductive volume or body 12' and the 
innermost surface of the annular support 22, i.e. the innermost surface of the agar jelly 

10 filling the compartments 32, to 32 16 . The sponge 36 is pliable and conforms to the shape 
of the conductive volume surface. It connects each of the electrodes to R 16 and to 
S 16 electrically to the surface of conductive volume 12'. The conductivity of the 
interfacing medium 36 is comparable to the conductivity of the imaging medium 12'. 

The mounting of the electrodes Sj to S 16 in compartments 32, to 32 16 is the same 

15 in each case, so only the mounting of electrodes S 4 and R 4 described - as an example. 
Thus, the stimulation electrode is mounted upon the inner face of band 30 and connected 
through the band to terminal 38 4 . The recording electrode R 4 is suspended in the middle 
of the compartment 32 4 by a plastics ligature 40 4 , extending from the annular flange 26 
and a conductive filament 42 4 extending from the other annular flange 24. The filament 

20 42 4 connects to terminal 44 4 on the outer surface of flange 24. The terminals 38! to 38 16 
and 44j to 44 4 are not shown in Figures 2 and 3 for clarity. 

The terminals 38 and 44 are connected to the electronic conditioning circuit 16'. 
For convenience, the individual connections are not shown in Figures 2 and 3. The 
circuit 16' permits selection of individual pairs of the outer electrodes Sj - S 16 for 

25 stimulation by a suitable current of, say, 0.1 - 1mA at 10 - 100kHz, and selection of 
pairs of the inner electrodes R x - R 16 for measurement of potential differences across 
them. Each of the electrodes R x to R 16 and S! to S 16 comprises a disc mounted with its 
plane normal to a radius of the annular support, each pair of electrodes in a cell being 
aligned radially, i.e. on a line extending normal to the surface of the body 12'. Each 

30 electrode is about 6 mm in diameter. The electrodes R t to R 16 and Sj to S 16 in their 
respective compartments are spaced apart radially by a distance sufficient to reduce/avoid 
shielding effects. Other sizes of electrodes, for example from 2 mm to 10 mm diameter, 
are available. If they were used, the spacing would be adjusted appropriately. 
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In view of the alternating "high" conductivity of compartments 32 x to 32 16 and 
"low" conductivity of separators to 28 16 , the annular support can be considered to 
have anisotropic conductivity, i.e. high conductivity for currents flowing radially from 
an electrode S towards the volume 12' and low conductivity for currents flowing 
5 circumferentially from one compartment towards its neighbour. 

The manner in which this arrangement of electrode geometry and media favours 
the enhancement of the sensitivity of the apparatus will now be described. 

Figure 5 illustrates the average sensitivity, in the radial direction, obtained using 
a finite element model, for an electrical impedance tomography apparatus using the 
10 annular electrode array of Figure 2, and that for an equivalent conventional electrical 
impedance tomography apparatus using the electrode array of Figure 1 . In both cases 
the conductive volume 12' has a radius of 10 cm, with a uniform conductivity. For 
proof of concept, the annulus in Figure 2 extends from 10 cm to 15 cm and has a 
conductivity 10 times larger than the conductivity of conductive volume 12'. The 
15 medium separating the individual compartments, i.e. of separators 28! to 28i 2 , is set to 
a conductivity value of 1% of the conductivity of conductivity volume 12'. 

The average sensitivity curve obtained for the conventional electrode array of 
Figure 1 uses configurations of stimulator and recording electrode pairs as disclosed by 
Brown et ah in U.S. patent No. 4,617,939. As mentioned earlier, in this approach, a 
20 pair of neighbouring electrodes is used to stimulate while the electrical potential is 
measured with the remaining pairs of electrodes. The process is repeated until all the 
independent combinations are exhausted. The electrodes used for stimulation, however, 
are not used also for measurement. 

The average sensitivity curve for the electrode configuration shown in Figure 2 
25 is obtained using pairs of electrodes in the outer group S r S 16 to apply stimulation current 
and the corresponding pairs of electrode in the inner group RrRi 6 to measure the 
resulting electrical potential. 

The resolution of the reconstructed image depends upon the sensitivity of the 
apparatus to changes of the electric conductivity in the conductive volume 12'. The 
30 sensitivities depicted in Figure 5 were obtained using the sensitivity expression described 
by D.B. Geselowitz in "An Application of Electrocardiographic Lead Theory to 
Impedance Plethysmography", IEEE Trans. Biomed. Eng., Vol. BME-18, pp. 38-41, 
Jan. 1971, and by J. Lehr in "A Vector Derivation Useful in Impedance 
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Plethysmographic field calculation", IEEE Trans. Biomed. Eng., Vol. BME-19, pp. 
156-157, Mar. 1972. 

The theoretical basis for the reconstruction algorithm employed by embodiments 
of the invention will be described with reference to Figure 6. When the conductivity 
5 distribution changes from a(x,y,z) to <r(x,y,z) + A<x(x,y,z), the change in the transfer 
impedance AZ for the pair of stimulating electrodes (A,B) and recording electrodes 
(C,D) can be given as: 



AZ = - 



V<|)(a) Vt|f(q + A a) 



dv 



(1) 



10 where: 

4> is the potential distribution when the current I* is applied to the 

electrode pair (A,B). 
^ is the potential distribution when the current I* is applied to the 

electrode pair (C,D). 

15 a(x,y,z) is the conductivity distribution when ^ is established for I* 

cr(x,y,z,) + Aa(x,y,z) is the conductivity distribution when ^ is 

established for I*. 

The term (cr+Acr) * n expression (1) can be expanded with respect to 

A<r. Assuming small Acr's, the higher order terms can be neglected and expression (1) 
20 is expressed as follows: 



AZ= -UaYm.Yl^dv 
J /„ /, 

(2) 

Since the term AZ, the difference between the measured transfer impedance and 
the assumed transfer impedance, is known, the difference in the conductivity distribution 
25 Acr can be obtained. Dividing the region of interest into small domains over which the 
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conductivity is assumed constant enables expression (2) to be discretized and expressed 
as the following matrix equation: 



M 

5 where: 

AZj is the difference in transfer impedance for the ith combination of a pair 

of stimulating electrodes and a pair of recording electrodes; 
Aoj is the unknown conductivity difference for the region element j; 
Sy is the sensitivity coefficient term for the stimulating/recording electrode 
10 pair combination i, at the element j; and 

m is the number of domains into which the region has been divided. 
The term Sy is given as follows: 



S„~\ V *' V * d Vj (4) 



15 An examination of expression (4) indicates that the sensitivity term can be 

enhanced by an apparatus that could maximize the integrand. Sensitivity enhancement 
is achieved by selecting configurations of stimulating and recording electrode pairs that 
favour the scalar product in the numerator and using an annulus medium with an 
anisotropic conductivity that limits the current intensities at the denominator. 

20 

Electrode configurations. The annular electrode array of Figure 2 allows electrode 
configurations that can maximize scalar products of stimulating and recording electric 
field vectors more so than could be obtained with the electrodes of Figure 1. It permits 
more combinations of electrodes to be made, with a greater number of different 
25 geometries. 

It should be appreciated that the stimulating electrodes need not be diametrically 
opposite each other. Referring to Figure 2, using pairs of electrodes in the outer group 
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Sj to S 16 to stimulate the body and pairs of electrodes in the inner group R t to R 16 to 
measure the resulting electrical potentials tends to maximize the scalar product in the 
numerator of expression (4). 

5 Annulus conductivity. The sensitivity of the apparatus to a change in the electrical 
conductivity at a position (x,y) in the conductive volume is an inverse function of the 
product of the total electric current intensities that would be required to produce the 
electric field intensity of the stimulating and recording electrode pairs at position (x,y). 
Disposing the electrodes in sectors whose boundaries are poorly conductive, produces 
10 an annulus medium with an anisotropic conductivity, vis. conductive in the radial 
direction and poorly conductive in the direction. This arrangement forces the current to 
flow in the medium of interest, i.e. the conductive volume, as opposed to the extended 
medium. It favours a larger ratio of the electric field to current intensity, in the medium 
of interest. 

15 The sensitivity function used to determine the average sensitivities represented in 

Figure 5 was not optimized for the various parameters that define the annulus. Key 
parameters are expected to be the following: 

A. Annulus dimension: The size of the annulus used for this evaluation has an outer 
20 to inner radius ratio of 1.5. Evaluations have shown that as the ratio increases, and the 

electrodes are set further away from the imaging volume, the sensitivity decreases. 
Conversely, as the ratio decreases, and the electrodes get closer to the imaging medium 
and the sensitivity increases. Although this is true for a discrete point electrode model, 
the optimal practical ratio will be determined by the size of the electrodes and the 
25 minimum separation required to avoid mutual shielding. 

B. Annulus electrical conductivity: Enhanced sensitivity is obtained when the 
annulus is configured with an anisotropic electrical conductivity distribution, that is, a 
large ratio of the radial to the azimuthal electrical conductivity. No significant 

30 improvement was gained by decreasing the azimuthal electrical conductivity below 1 % 
of the average imaging medium electric conductivity. 
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C. Electrodes: A minimum of 32 electrodes were considered, arranged in two 
concentric rings. For this evaluation, each wedge-shaped compartment held two 
electrodes, aligned along the diameter. One electrode was positioned at the edge of the 
annulus and the other was positioned at the centre of the annulus. Although this was 
5 effective for the model, in practice, of course, the other electrode cannot be any closer 
to the centre than the inner surface of the annulus. It is expected that standard 
commercial AgCl disk shaped electrodes would be used for this application. 

Measurement procedure. As stated earlier, it is desirable to maximize the scalar 

10 product of the electric field vectors for the stimulating pair of electrodes and the 
"recording" electric field vectors for the recording pair of electrodes. Consequently, for 
a given pair of stimulating electrodes, it is sufficient to measure the potential difference 
between that pair or pairs of electrodes which are closest spatially to the stimulating 
electrodes and not to measure the potential differences between other pairs of electrodes. 

15 This measurement procedure will now be described with reference to the flowchart of 
Figure 7. In order to illustrate the stimulation/measurement strategy explicitly, an 
example will be given using the geometry depicted in Figure 2. The outer electrode set 
(SrS 16 ) are used as stimulating electrodes and the inner electrode set (R r R 16 ) as 
recording electrodes. It should be appreciated that the outer set could be replaced with 

20 a set on the inner surface, as per the comments above. 

Following positioning of the electrodes S r S 16 and RrR 16 upon the conductive 
volume to be imaged, an initial conductivity distribution appropriate to the volume would 
be assumed (step 7.1). In the absence of other information from, say, previous 
measurements of the same or similar volumes, a uniform mean conductivity distribution 

25 would be assumed. For example, for landmine detection, this might be a typical 
conductivity distribution of the particular type of soil whereas, for human body imaging, 
it might be a typical conductivity distribution for that part of the human body. 

In step 7.2, a series of measurements are made, passing a known current through 
each pair of stimulation electrodes S Y to S 16 in turn and recording the potential difference 

30 across the "closest" pair of electrodes, i.e. having the recording electric field vectors 
most closely aligned with those of the stimulating pair of electrodes. For example, when 
stimulating Sj and S 2 , the potential difference between Rj and R 2 only would be 
measured. 
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The stimulation/recording pairings proposed here are: 



Stimulation Pair 


Recording Pair 


Pair Gap 


Si 5S2 


Ri ,R2 


1 




R2,R3 


1 


J 9 *t 


jv*%Ra 


1 








Si, s 3 


1*1 ^3 


2 


S25S4 




2 








Si,S 4 


Ri,R 4 


3 









Table 1 Stimulate/Measure Configuration 



15 This pattern yields 120 independent combinations for 16 compound electrodes 

(stimulate-record pairs). It should be noted that the stimulating electrodes and the 
recording electrodes are always identically numbered, and hence closely located. Other 
stimulation/recording electrode configurations which share this feature will have some 
benefit. An example of such an alternative, possible with only a single set of electrodes 

20 as shown in Figure 1 (denoted E^E^ ... E 16 ) is: 



25 



Stimulation Pair 


Recording Pair 


Pair Gap i 


Ei,E3 


E2>E 4 


2 


Er2,E 4 


E3,E 5 


2 


E3,E 5 


E 4 ,E 6 


2 








Ei,E 4 


E2,E 5 


3 


E2>E 5 


E3,E 6 


3 









30 Table 2 Alternative Stimulate/Record Configuration 



This pattern has a reduced number of independent combinations, since it starts 
with a stimulating electrode gap of 2, for example, but will show some of the benefit of 
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the proposed configuration, and is shown as a simple example of an alternative 
manifestation of the stimulate/record configuration described herein. 

Referring again to the flowchart of Figure 7, once the measurements have been 
made, in step 7.3, the region of interest, i.e. the region in which current might flow, is 
5 divided into a mesh of elements. In Figure 2, the region of interest is the conductive 
volume 12', the sponge 36 and the annular array of compartments 32 r 32 16 . Each 
element has dimensions equal to approximately 10 per cent of the spacing between each 
pair of the recording electrodes R x to R 16 . The mesh could be triangular or rectangular, 
as desired, and the region could be two-dimensional or three-dimensional. 
10 In step 7.4, the sensitivity matrix is calculated for the mesh using the electrode 

geometry, initial conductivity distribution, and applied currents, according to equation 
[4], for the region of interest. 

In step 7.5, an estimate is made of the noise present in the potential difference 
measurements divided by the applied stimulation current. This scaled noise value is used 
15 to determine a parameter /jl which is the ratio of the scaled noise value to the root mean 
square (RMS) of the initial conductivity distribution. 

In step 7.6, the pseudo-inverse S # y of the sensitivity matrix is calculated 
according to equations (5) and (6). 

In step 7.7, the differential conductivity distribution Aa is calculated according 
20 to equation (7). 

Preferred electrode configuration and compound electrodes. 

Compound electrodes are so named because they comprise two parts, one part 
being used for stimulation and the other for recording. It has been found that compound 

25 electrodes, in conjunction with the novel measurement procedure described above, 
provide significantly better performance in physically realistic (noisy) situations. 
Analysis has shown that this improvement in performance is due both to sensitivity 
amplitude enhancements (due to alignment of the electric field stimulating and recording 
electric field vectors) and a better representation (more suitable for reconstruction) of the 

30 set of measurements possible by the proposed recording configuration. 

A potential problem associated with compound electrodes is the possibility of a 
partial electrical short circuit between the stimulating and measuring electrode, due to 
the close proximity of these electrodes. Simulations have indicated that satisfaction will 
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be obtained with spacing of the component electrodes of a compound electrode such as 
to render the voltage associated with this extra (surface) current path approximately 0.1% 
of that associated with the conventional path (through the region of interest). For a 
human body this requirement translates into a resistance of (for skin surface alone) 
5 approximately 5MQ. An example of a suitable compound electrode is shown in Figures 
8 and 9. The compound electrode comprises an annular outer electrode 46 surrounding 
a circular inner electrode 48 with an annular slot 50 between them. The electrodes 46 
and 48 are bonded to an insulating planar support 52 and their exposed surfaces are 
coated with layers of conductive gel 54 and 56, respectively. An "O" ring 58 of 

10 insulating material is located in the annular slot 50. The diameter of the "O" ring 58 is 
approximately equal to the depth of the annular slot 50 so that, when the compound 
electrode is applied to the surface of the volume under investigation, the "O" ring 58 will 
prevent gel 54/56 from exuding across the annular slot 50 and causing, effectively, a 
short-circuit between the inner and outer electrodes. In use, a plurality of such 

15 compound electrodes would be distributed about the surface of the body being 
investigated. Current could be applied to the respective outer electrodes of a pair of 
compound electrodes and the potential difference recorded between the inner electrodes 
of the same pair. Conversely, the inner electrodes could be stimulated and the potential 
difference recorded between the outer electrodes. 

20 It should be appreciated that other configurations of compound electrode could 

be used. For example, either or both of the electrodes could be rectangular. It is also 
envisaged that each compound electrode might comprise an array or grid of stimulation 
electrodes interspersed with an array or grid of recording electrodes. 

The annular electrode array of Figures 2 and 3, whether with the individual 

25 electrodes or with the compound electrodes, is well suited for medical diagnostics 
applications applied to a patient's head, limb or torso. Nevertheless, the enhanced 
sensitivity approach illustrated with the annular geometry can also be applied to other 
electrode configurations. Thus, Figure 10 illustrates a two-dimensional array of 
electrode cells which is shown as rectangular but could have other shapes. This 

30 electrode array is envisaged for ground imaging applications, such as landmine detection. 
The array comprises juxtaposed individual electrode cell units 60. As shown in Figure 
11, each cell comprises a pair of electrodes R, S immersed in a conductive material, for 
example a jelly formed from agar and saline solution, as before. 
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The cell units 60 lie on a bed 62 (Figure 10) of conductive material that provides 
the electrical interface with the ground. This bed 62 would have insulating partitions 
(not shown) to maintain anisotropic conductivity. 

Figure 12 illustrates another possible embodiment, where a flexible vest 64 
5 comprises similar individual electrode cells 66 forming a cylindrical array which, as 
shown, can be wrapped around a human torso and used for three-dimensional imaging. 
It will be appreciated that other geometries are also feasible, such as a special cap for 
head imaging applications or a special brassiere for breast imaging applications. 

It should be appreciated that, if the radially-spaced electrodes S/R of Figures 2-4 
10 and 10-12 were replaced by compound electrodes, the agar-filled compartments or cells 
would not be required, since the compound electrodes could lie flat upon the surface. 
However, it is envisaged that, to facilitate calculations, the compound electrodes would 
still be in a rigid support, if necessary, with an intervening sponge 36, or the like. 

15 Reconstruction technique 

The proposed reconstruction technique is based on the construction of a 
probabilistically regulated, pseudo inverse of the sensitivity matrix introduced above and 
mentioned in step 7.4. The construction of this pseudo inverse is very straightforward, 
and follows the following steps :- 
20 Calculate the singular value decomposition of the sensitivity matrix. 

S — U A V T (5) 

Replace each singular value \ in A = diag(X,) with x* =;t. / (X] + \x 2 ), 

producing A # = diag{k]) . The parameter fi is the ratio of the noise in the measured 

potential differences to the root mean square (RMS) of the conductivity distribution. 
25 Thus, in this sense it is the inverse of the SNR. 
Form the pseudo inverse 

S» = VA« U T (6) 
which can be used to compute conductivity perturbations from potential difference 
recordings, i.e. reconstruct the EIT image, since 

30 
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Ao i = H S ij AZ j> < 7 > 

where n is the number of independent measurements taken. 

This approach has been found to respect the inherent low-pass filtering of the EIT 
process, and offers a rapid and robust option for EIT image reconstruction. 
5 Figures 13 A to 13J illustrate a comparison between the novel measurement 

procedure set out in Table 1, using compound electrodes as per Figure 8, and the 
measurement procedure disclosed by Brown et ah in U.S. Patent No. 4,617,939, but in 
both cases using the reconstruction technique described above. Figures 13A to 13J 
depict the reconstruction of a perturbation covering approximately 10% of the diameter 

10 of a circular region of interest at three different locations. As shown in Figures 13A, 
13D and 13G, with zero noise, both Brown et al. 's procedure and the procedure of the 
present invention reproduce the perturbation very accurately. Thus, Figure 13A shows 
a maximum a posteriori (MAP) estimate of a perturbation near the circumference of the 
region. Figure 13B shows the (MAP) estimate of the differential conductivity 

15 distribution A<r for a circular array of 16 compound electrodes of the kind shown in 
Figure 8. Figure 13C shows the corresponding MAP estimate for the conventional 
electrode configuration of Figure 1 using Brown et al 's measurement technique. 

Figures 13D to 13F correspond to Figures 13 A to 13C for a perturbation which 
is closer to the middle of the region and Figures 13G to 13J are the corresponding 

20 reconstructions for a perturbation at the centre of the region. 

Although embodiments of the invention have been described and illustrated in 
detail, it is to be clearly understood that the same are by way of illustration and example 
only and not to be taken by way of the limitation. The spirit and scope of the present 
invention embrace various modifications and substitutions. For example, the electrodes 

25 could be in a fixed array and a separate current generator provided for each pair 
stimulated. Alternatively, a single current generator could be used to stimulate each pair 
of electrodes in turn, using a commutation switching circuit, enabling a more expensive 
current generator to be used. 

It is also envisaged that only two stimulating electrodes and two recording 

30 electrodes could be used. The various locations would be predefined and the electrodes 
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repositioned at selected locations in turn until all locations had been stimulated as 
required. 

It is envisaged that, instead of providing a personal computer with a DSP card 
for performing the reconstruction, a suitable custom integrated circuit, preprogrammed 
5 with the reconstruction algorithm software, could be combined with the electronic 
conditioning circuitry in an interface unit which would connect to a conventional personal 
computer and/or even directly to a display device. 

INDUSTRIAL APPLICABILITY 
10 The compound electrode arrangement may be easier to manufacture. The 

measurement procedure facilitates simpler implementations for a given sensitivity. 



WO 98/23204 



PCT/CA97/00876 



20 



CLAIMS: 



1. A method of determining electrical impedance tomography of a conductive 
volume (12) comprising the steps of selecting, in turn, a plurality of pairs of locations 
5 (S) upon a surface of the volume (12) and supplying a stimulation current to the surface 
by way of the selected pair of locations, and recording, for each selected pair, a resulting 
potential difference between at least one pair (R) of the remaining locations, and 
processing the recorded potential differences to plot electrical conductivity distribution 
within said volume characterized in that , for a particular pair of locations stimulated, the 
10 potential difference is recorded for a number of pairs of said locations that is less than 
the total of said plurality of locations and also comprises that pair of locations which, if 
stimulated, would produce an electric field with vectors most closely aligned with the 
corresponding vectors of the electric field produced by stimulation of said particular pair 
of locations. 



2. A method of determining electrical impedance tomography according to claim 1, 
using a plurality of electrodes (Sj-S^, RrRi 6 ) disposed upon or adjacent said each surface 
at a respective one of said plurality of locations, characterized by the steps of; 



15 



20 



(i) 



applying a stimulation current to said volume by way of a first and a 
second of said plurality of electrodes at corresponding first and second 
locations; 



25 



(ii) 



simultaneously recording the resulting potential difference between that 
pair of the remaining electrodes at said pair of locations which, if 
stimulated, would produce electric field vectors most closely aligned with 
the electric field vectors produced by stimulation of said particular pair of 
locations; 



(iii) 



applying a stimulation current to said body by way of two electrodes other 
than those stimulated in a previous step; 



30 



(iv) 



simultaneously recording the resulting potential difference between that 
pair of the remaining electrodes at said locations which, if stimulated, 
would produce an electric field with vectors most closely aligned with the 
corresponding vectors of the electric field produced by the stimulation 
current applied in step (iii); 
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(v) repeating steps (iii) and (iv) until all of said plurality of electrodes have 

been stimulated and the corresponding potential differences recorded. 

5 3. A method according to claim 2, characterized in that stimulation current is 
applied to a said first and second of said plurality of electrodes with no intervening 
electrodes. 

4. A method according to claim 2, characterized in that said stimulation current is 
10 applied to a said first and second of said plurality of electrodes that are spaced apart with 

one or more other electrodes therebetween. 

5. A method according to claim 4, characterized in that , each performance of steps 
(iii) and (iv) applies said stimulation current to a pair of electrodes that are spaced apart 

15 with a different number of intervening electrodes therebetween. 

6. A method according to claim 2, characterized in that , in said steps (i) and (ii), 
said stimulation current is applied to a said first and second of said electrodes that have 
no electrodes therebetween, and said potential difference recorded, and, in said steps (iii) 

20 and (iv), a stimulation current is applied to a said two electrodes that have one or more 
electrodes therebetween, and the resulting potential difference recorded. 

7. A method according to claim 2, characterized in that the stimulating currents are 
applied to selected pairs of a first group of the plurality of electrode means and the 

25 resulting potential differences are recorded at selected pairs of a second group of the 
plurality of electrodes. 

8. A method according to claim 7, characterized in that the electrodes of one of the 
groups are arranged in a ring closer to the surface than the electrodes of the other of the 

30 groups, each electrode of the first group being substantially aligned with a corresponding 
electrode of the second group in a direction extending substantially normal to said 
surface and said potential difference is recorded between that pair of said electrodes 
which are aligned with the pair of electrodes to which the stimulation current is applied. 
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9. A method according to claim 7, each electrode means comprising a stimulation 
electrode and a recording electrode closely located to each other and both the same 
distance from said surface, characterized in that the step of applying stimulation current 
applies the stimulation current to the stimulation electrodes of two electrode means, 

5 respectively, and the corresponding potential difference is measured between the 
recording electrodes of the same two electrode means. 

10. An electrode arrangement for an electrical impedance tomography system 
comprising a first group of electrodes (R r R 16 ) and a second group of electrodes (S r S 16 ) 

10 and support means (24,26,28,30) for supporting the electrodes adjacent a surface of a 
volume the electrical impedance of which is to be mapped by the system, characterized 
in that the first group of electrodes is closer than the second group of electrodes to a part 
of the support means which, in use, will be adjacent said surface. 

15 11. An electrode arrangement according to claim 10, characterized in that each 
electrode of the first group (R r R 16 ) is substantially aligned with an adjacent electrode of 
the second group (S r S 16 ) in a direction extending substantially normal, in use, to said 
surface. 

20 12. An electrode arrangement according to claim 10 or 11, characterized in that the 
support means comprises an annular structure, the first group of electrodes and the 
second group of electrodes being arranged in an inner ring and an outer ring, 
respectively, about a central part of the annular structure. 

25 13. An electrode arrangement according to claim 12, characterized in that the annular 
structure is circular and the spacing between the inner ring of electrodes and the outer 
ring of electrodes is about one quarter of the diameter of the inner ring. 

14. An electrode arrangement according to claim 10, characterized in that the support 
30 means comprises a plurality of compartments (32 r 32 16 ), each comprising a conductive 
medium for contacting the surface of a volume to be imaged, and lower conductivity 
partitions (28 r 28 16 ) between adjacent compartments, each compartment housing a pair 
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of electrodes comprising one electrode from said first group and one electrode from said 
second group. 

15. An electrode arrangement according to claim 14, characterized in that the 
5 electrodes are supported in a planar array (60,62), and each electrode of said first group 

is spaced from an associated electrode of said second group in a direction substantially 
perpendicular to the plane of the array. 

16. An electrode arrangement according to any one of claims 10 to 15, characterized 
10 in that the support means is rigid and a conductive interface medium (36) is disposed for 

juxtaposing between an inner surface of said support means and the surface of the 
volume (12) to be measured, the interface medium comprising a material that is 
conformable to contour variations of said surface of the volume. 

15 17. An electrode arrangement according to claim 16, characterized in that the 
interface medium material has a conductivity comparable to that of the volume to be 
measured. 

18. An electrode arrangement according to any one of claims 10 to 17, characterized 
20 in that the support means has anisotropic conductivity with its conductivity in said 

direction normal to said surface, in use, being significantly higher than its conductivity 
in directions transverse thereto. 

19. An electrode arrangement according to claim 18, characterized in that the support 
25 means comprises a plurality of higher conductivity sectors alternating with, and defined 

by, a plurality of lower conductivity regions. 

20. An electrode arrangement for an electrical impedance tomography system 
comprising a plurality of electrodes (S r S 16 , RrRi 6 ) in an array mounted in a support 

30 (24,26,28,30) for supporting the electrodes adjacent a surface of a volume (12') to be 
imaged, characterized in that the support medium has anisotropic conductivity, its 
conductivity in a direction that, in use, is normal to the surface, being significantly 
greater than its conductivity in a transverse direction. 
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21 . An electrode arrangement according to claim 20, characterized in that the support 
comprises a plurality of compartments (32 r 32 16 ), each separated from neighbouring 
compartments by an insulating partition (28 r 28 16 ), each compartment comprising a 
material having a conductivity comparable to the conductivity of the volume the 

5 conductivity distribution of which is to be measured and containing at least one of the 
electrodes. 

22. An electrode arrangement for an electrical impedance tomography system 
characterized by an array of electrodes fixed spatially relative to each other in a rigid 

10 support, and an interface medium for interfacing the rigid support to a surface of a 
volume the conductivity distribution of which is to be measured, the interface medium 
being conductive so as to connect the electrodes electrically to the surface and pliable so 
as to conform to variations in relief of the surface. 

15 23. An electrical impedance tomograph comprising an electrode arrangement 
according to any one of claims 10 to 22, further characterized by means (16,18,20) for 
stimulating, sequentially, pairs of electrodes and, for each pair, measuring corresponding 
potential differences between pairs of the remaining electrodes, and processing the 
measured potentials to map electrical conductivity distribution of a volume. 
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